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^ ■ Abstract 

We have calculated the leading Yukawa corrections to the chargino, neutralino 
t^^ \ and gluino pole masses in the DR scheme in the Minimal Supersymmetric Standard 

^D ■ Model (MSSM) with the full set of complex parameters. We have performed a 

numerical analysis for a particular point in the parameter space and found typical 
corrections of a few tenths of a percent thus exceeding the experimental resolution 



^ ■ as expected at the ILC. We provide a computer program which calculates two-loop 

^ ■ pole masses for SUSY fermions with complex parameters up to 0{aas-,Y'^,asY'^). 



1 Introduction 

In the Minimal Supersymmetric Standard Model (MSSM), each Standard Model particle 
is part of a SUSY multiplet containing its superpartners of opposite statistics and a spin 
differing by 1/2 unit. Moreover, due to anomaly cancellation the SUSY Higgs-sector 
consists of two chiral multiplets whose fermionic degrees of freedom are called higgsinos. 
On the other hand, the vector supermultiplets contain the standard model gauge bosons 
and their spin-1/2 superpartners called gauginos. 

Organizing these particles in mass eigenstates one finds four mixing neutral Majorana 
fermions composed of the superpartners of the photon, the Z^-boson, and the neutral 
Higgs bosons -ff? 2- Then there are two charginos Xi and X2 ; which are the fermion mass 
eigenstates of the partners of the W^ and the charged Higgs bosons Hi2- Finally, the 
gluino is the superpartner of the SU{3)c vector-boson, the gluon. 

The next generation of future high-energy physics experiments at Tevatron, LHC 
and a future e~^e~ linear collider (ILC) will hopefully discover some of these particles 
if supersymmetry (SUSY) is realized at low energies. The physics interplay of these 
experiments has been studied in [1] and it turned out that particularly a linear collider 
[H [2], [3] will allow to test the underlying SUSY model with great accuracy. In fact, the 
accuracies of the masses of the lighter SUSY fermions are expected to be in the permille 
region which makes the inclusion of higher order loop-corrections indispensible. 

Within the real MSSM important results on quark self-energies were obtained in [1]- 
[6]. In [Z]-[in| the gluino pole mass was calculated to two-loop order 0{a'g). Moreover, 
the scalar MSSM Higgs-sector has been studied in detail, also in the full complex model 
[TT]-[T3] and even some three- loop results have been included recently [TB] . 

In previous works [TTJ [18] we studied the SUSYQCD two-loop corrections to chargino 
and neutralino pole masses and found that these effects have to be taken into account 
when matching DR -input with experimental data. However, the reference scenario used 
for the numerical analysis was the benchmark point SPSla' [12], a scenario where the 
lighter mixing partners have a dominating gaugino component. But once the higgsino 
component becomes more important one can expect that also leading two-loop Yukawa 
corrections are of the same size. 

In this paper we therefore study these leading two-loop Yukawa corrections to the 
pole masses of charginos, neutralinos and the gluino. Since the one-loop correction to the 
gluino pole mass is of order 0{as) the leading two-loop Yukawa correction is of order 
0{asY^). On the contrary, leading Yukawa corrections for charginos and neutralinos are 
of order 0(Y^). We neglect the Yukawa coupling to leptons. The reference point for the 
numerical study is chosen such that the lighter mixing partners have a dominant higgsino 
component and their masses are roughly of the size of their SPSla'-values so that we can 
compare the magnitude of the corrections to the expected experimental accuracy at this 
point. 

We conclude that the leading two-loop Yukawa corrections typically exceed experimen- 
tal resolution for the case of neutralinos. For charginos we find smaller corrections com- 
parable to their two-loop SUSYQCD corrections [18] which should be taken into account 



when extracting DR parameters from experiment. The leading two- loop Yukawa correc- 
tions to the gluino are comparably small but extending the known two-loop SUSYQCD 



corrections [9l [10] to include the gluino phase for the numerical evaluation gives important 
results. 

For the calculation of the self-energy amplitudes we employ semi-automatic Math- 
EMATICA tools [201 EI] and cross-check the results with the generic analytic formulae 
for two-loop corrections to fermion pole masses in a general renormalizable field-theory 
derived in [9]. 

In principle, some of the complex phases of the full complex MSSM can be rotated 
away by rigid field transformations, e.g. the phase of the gluino mass parameter. However, 
we do not pose restrictions on these phases for convenience. Instead we generalize the 
explicit results on the SUSYQCD corrections to SUSY fermions derived in [9] to the case 
of a complex gluino mass parameter. For the above reason this is not a generalization of 
the physics. 

All of these new results are implemented in the C-program POLXINO 1.1 [22] • It 
currently is the only publicly available program for the calculation of even the one-loop 
chargino and neutralino pole masses with complex parameters and described in |18] . 

We list all the analytic results in the appendix. 

2 Diagrammatics and calculation 

The masses of a set of mixing fermions are defined through the complex poles of the full 
renormalized propagator of the mixing system as gauge-invariant and renormalization 
scale invariant quantities [23]- [31]. The tree- level values are the eigenvalues of the tree- 



level DR mass-matrix m. 

This complex pole masses can be inferred from the zeros of the determinant of the 
renormalized effective action T^'^^p'^) in terms of the renormalized self-energies S 

-^f(2)(p) = G(2)-l(p) 

[a-pil + tyip^]) -mt + S^/(p2) )■ 

(1) 

Here, we use the usual definition for the Lorentz-generators a'^ = (l,cr*) and a^ = 
(1,-0"*). The diagonalized tree-level mass-matrix m could be chosen real and semi- 
positive definite by a rigid reparametrization of the fields. However, in the case of a 
complex gluino mass parameter we find it more natural to keep its complex phase in the 
bilinear Lagrangian and therefore we keep the complex conjugates in Eq. ([T]). 

The pole mass condition in the Weyl representation ([1]) becomes particularly simple 
in the rest frame where kinematics is trivial 

= det(s - (m - S^(.)) ■ (1 + ti{s)r' ■ (mt - t^{s)) ■ (1 + tf{s)r'). (2) 

Note that the objects S^'-? are finite matrix valued functions of the external momentum 
invariant s = p^ . The iterative solution to Eq. ([2]) is [9] 
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All self-energies have to be evaluated at an external momentum invariant s = |mjp with 
an infinitesimal positive imaginary part |32j . 

According to the DR -scheme the UV-divergencies of the Feynman-integrals are reg- 
ulated dimensionally hj d = 4 — e, and the unphysical scalar mass parameter rrig for the 
evanescent fields is absorbed according to [5B] . 

The resulting two-loop integrals have a rather complex tensor structure which has 
already been worked out in detail in [S], implemented in a C-program [TU] and tested 
thoroughly [171 [IB]. In this work we therefore use the notation of [9] for the basis integrals. 

2.1 One- loop level 

In Fig.[T}l3]we show the one-loop diagrams for the gluino, the neutralinos and the charginos, 
respectively. Note, that diagrams with charge conjugated inner particles are not shown 
explicitly. We checked our analytic one-loop calculation against previous work [M| 135] in 
the on-shell scheme and found agreement. 

In the case of the gluino the one-loop correction is proportional to g'^, thus any squared 
one-loop self-energies cannot contribute to leading two-loop Yukawa corrections 0{g'^Y'^). 
The result is 

U^W = ^^2g2 {c^\mf (3E|m~|2 - 5) - Bpsiq, Qs) + 26^^(9, q,)mgRe{m~Ll* Rl)) 

(4) 

which slightly generalizes Eq.(5.2) of [9] for a complex gluino mass parameter. 

The 0{a) one-loop mass shifts for charginos and neutralinos are given in Eq.(5.15) and 
Eq.(5.18) of [9] and are not repeated here. However, for a fixed order calculation 0{a, Y"^) 
we need the 0(Y^) one-loop self-energies because their squares appear in Eq. ((HI). 

The results are 

^m ij = (^m ) ^j = ~ J^^qxJ Qt^x^ qQ s^F S il ^ Q s) 

^l- = ^— ¥*-_-¥-._- BFsiq,Qs)/s 

16tt^ i^iQs iXjQ*s ■f'-sv'i' ^«^/ 



167r2 Xj qQs Xi iQs ^^^'^^'^^J/ 
^m ij = (^m ) ij = ~Jq:^ \^gx°qt^x°qqs + ^qx°q*s^x°qqs J ^Fsi^^ Qs) 

sf .■ = (Ef )% = -^ (>~o,~ r^o,-~ + i;~o~.i;^o~,) bMq, Qs)/s (5) 

where summation over the squark index s and q = {u, d, c, s, t, b) is understood imphcitly. 
Q denotes the iso-spin partner of q and is not summed over independently. The Yukawa 
couphngs are given by 
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where we employ the conventions of [36] for the mixing parameters. Furthermore, 

Yu = —F^ and Y^ = —;= . (7) 

In Eq. ([6]) and Eq. ([7]) the symbols u and d denote the up-type and down-type quark, 
respectively. The generation index is suppressed. 

2.2 Two-loop results 

In Fig.Hland Fig. Owe show the leading two-loop self-energy diagrams contributing 0{Y^) 
to the pole mass. For simplicity we put the external neutralinos and the gluino in Fig. H] 
into a generic Majorana field 77. Without loss of generality neutralino and chargino masses 
are taken to be real, the gluino mass however is allowed to have a complex phase thus 
avoiding the introduction of a one-dimensional rotation matrix for this particle. 

The renormalization scheme adopted is the familiar DR -scheme which regulates UV- 
divergencies dimensionally but introduces an unphysical scalar field for any gauge field 
in the theory in order to restore the counting of degrees of freedom in supersymmetry. 
These unphysical mass parameters can be absorbed in the sfermion mass parameters [33] , 
the resulting scheme is called DR'. 

The diagrams and the amplitudes were generated in FeynArts 3.2 [20] and simplified 
analytically using FeynCalc 4.0.2 [21]. The resulting tensor integrals were translated 
into the conventions of [9]. We then independently derived the same results from the 
generic formulae given in Eq. (4.4, 4.5) of [9]. 

The arguments of these self-energy tensor-functions are the squared masses of the 
inner particles and we abbreviate them by their respective symbol. The external mo- 
mentum invariant is the squared tree-level DR -input mass-parameter s = |mjjp with an 
infinitesimal positive imaginary part. 



For the numerical evaluation of the tensor integrals in our program Pobdno [22] we 
use the implementation of [10] and the two- loop self-energy library [32] . 

We split the rather lengthy analytic expression into five parts according to which 
particles occur in the loops, neutral (scalar and pseudoscalar) Higgs bosons, charged 
Higgs bosons, neutralinos, charginos or only quarks and squarks 

nP^ = (T^(%i?o + n~|^-+n~|~+n~|~o + n~|~.) (9) 

where x iii Eq. (Q can be either the neutralino or the chargino. The explicit form of 
the self-energy functions on the r.h.s. are given in the App. (]A.ip - (]A.3l) . Due to charge 
conservation there are quark isospin-partners in the chargino diagrams in the loop (Fig. [5]) 
denoted by {q,Q)- In order to get the leading Yukawa correction (9(y^) at the two-loop 
level we have to use the Yukawa part of the couplings given in Eq. (^ and Eq. 
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Figure 1: Gluino one-loop self-energy diagrams 
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Figure 2: Neutralino one-loop self-energy diagrams 
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Figure 3: Chargino one-loop self-energy diagrams 
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Figure 4: Leading Yukawa gluino and neutralino two-loop self-energy diagrams, t] = 
idyX^)- Q denotes the isospin partner of q = {u,d). cfp labels the neutral Higgs scalars 
{h'^ , H^ , G^ , A^) and (j)~ the charged ones {G~,H~). The direction of the charge flow is 
consistent with q = u. The arrows on quark and squark lines have to be reversed for q = d. 

3 Numerical results 



Our reference scenario for the numerical analysis is a higgsino scenario defined through 
Ml = 200 GeV, Ma = 400 GeV and /i = 130 GeV at Q = 1 TeV. All other parameters 
are taken from SPSla', in fact tan^ = 10, At = -532.38 GeV, A = -938.91 GeV, 
Mq3 = 470.91 GeV, M(j^ = 385.32 GeV and M^^ = 501.37 GeV, for further details see [19]. 
Evaluating the DR -spectrum, the one-loop corrections and our two-loop results Eq. ([HI 
[9]) as well as the respective two-loop SUSYQCD corrections [9] we find the following result 
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Figure 5: Leading Yukawa chargino two-loop self-energy diagrams. Q denotes the isospin 
partner of q = (u,d). cfP labels the neutral Higgs scalars {h'^ , H^ , G'^ , A^) and 0~ the 
charged ones {G~ ,H~). The direction of the charge flow is consistent with q = u. The 
arrows on quark and squark lines have to be reversed for q = d. 



for the real parts of the one- and two-loop masses : 



A.4 



106.54 -0.86 -O.OSGeV 
137.22 +1.92 +0.15GeV 
212.81 -5.66 + O.OSGeV 
417.87 -0.40 +0.04GeV 






121.54 +0.86 -0.05GeV 
417.75 -0.09 +0.02GeV 
572.33 +39.1 + 10.9 GeV 



(10) 



In this scenario the lighter mixing partners of charginos and neutralinos (x?, X2' XF) 
have a dominating higgsino component but their tree-level masses are roughly of the size of 
the SPSla' values. Therefore, we can compare the magnitude of the two-loop corrections 
to the expected experimental accuracy at SPSla'. 

In Fig. [6] we show the neutralino pole masses m^o and m^o and the two-loop mass shifts 
on the right as a function of tan/3. The dotted line contains the two-loop SUSYQCD 
corrections only whereas the solid line is the full Olag^,^^) result. Despite the obvious 
cancellation between SUSYQCD- and Yukawa corrections for x^ the corrections exceed 
the expected experimental accuracy of about ±0.05 GeV [Tj. 

For Xi in Fig- El we find corrections of a few hundred MeV for rather high values 
of tan/3 which should be included when matching DR-parameters to on-shell masses. 
However, corrections to Xi tend to stay below expected experimental accuracy. 

From the dependence of 6^'^^m^o and 5^^'m^o on At in Fig. [8] we learn that At severely 



effects these two-loop corrections and typical results exceed experimental resolution by an 
order of magnitude. On the contrary, S^'^^m-- in Fig. [9] as a function of At is again smaller 
than the expected experimental accuracy. In Fig. [TU] we show the pole masses and the 
corrections as a function of (pAt for the two light neutralinos. The two-loop corrections 
substantially depend on this phase. 

Note that for Xi ^iid Xi the Yukawa corrections tend to balance the SUSYQCD 
corrections. 

In Fig. [TT] and Fig. [T2] we investigate the dependence on (pMs of the lightest mass- 
eigenstates and the gluino. For the gluino it turns out that only the variation of the one- 
loop correction is numerically significant. Note also, that the two- loop Yukawa corrections 
to the gluino pole mass are very small. On the other hand, S^'^^m^o and 6^'^^m-- have a 
comparably strong dependence on 0^3 coming from the SUSYQCD corrections alone. 

4 Conclusions 

We have calculated leading two-loop Yukawa corrections to the pole masses of charginos, 
neutralinos and the gluino. The magnitude of these corrections typically exceeds exper- 
imental resolution for neutralinos. For charginos we find smaller corrections which still 
need to be taken into account when analyzing precision experiments. The leading two- 
loop Yukawa corrections to the gluino pole mass are rather small but we extended existing 
results on the SUSYQCD corrections to the case of a complex gluino phase giving rise to 
numerically relevant one- and two-loop effects. 

Finally, we included all results in the publicly available program POLXINO 1.1 
which has an convenient interface to the commonly used SUSY Les Houches accord 
for further numerical studies. 

A Analytic results 

In the following q is summed over {u,d,c,s,t,b) and Q denotes the iso-spin partner of 
q. The chargino index c, the neutralino indicex n as well as the squark-indices s, t, u are 
summed over independently. The generation index g appearing in Eq. (TT4l) is summed 
over (1, 2, 3). Finally, cfP is summed over (/i°, i7°, yl°, G°) and 0^ over {H~ ^ G~). |0°| is 
for the scalar {h^^H^) and -|-1 for the pseudo-scalar Higgs [A^fi^). The tensor-integral 
functions are listed in [9], their arguments denote the squared DR -masses of the particles. 

A.l Gluino 

For the gluino we get 

n^^'^ = (T^(%i?o + %|^-+n~|~+n~|~o + n~|~_) (ii) 

with 

ng|0o = g] (4mlYgg^oRe{m~X^Og^~.R'i*Ll)Mggj^g{q„ qt, q, q, 0°) 

9 



+ 4 Y,g^oRe{m^X^o~~,RtLr)MssFFs{qs, Qt, q, q, cfP) 
-Am,Y^,^.Re{\^,~~.{LfLl + {-)\^''\RTm)nssFFs(sis.quq.qA') 
+ 4 mgRe(m~A^Og^~, A*o~~* iifi^t ) V;p5S55(g, g„ g*, g„, <f) 

- 2 1 A<^o~~* pVir5S5s(g, gs, gs, gi, 0°) 

- 2 mqRe{m~X^a^o,~~,RfL'l)Ypggg{q, g„ g^, /) 
+ Re(A^o<^o.~~.)YF5S5(g, g„ g„ 0°) 

+ 4 (-)l<^Vjr^2_^oRe(i?^Lrmi) V5^p^5(g„ g, g, g, cfP) 

- 2 mJi;|^o(V5^^^ + 2{~t'\Y sF-WFs)isis. ?, ?, ?, /) 

+ Am,Y^^^oKe{m^RlLT){2NsF-FFs + (-)''^°'V5^^^5)(g., g, g, g, 0°) 

- 2 Y^g^oYsFFFsisis. g, g, g, 0°)) (12) 



n^i^- = ^.X4mgmQRe(m~rgQ^-A_^_,Q~.i?^Li^*)Mg5pp5(g„Qi,g,Q,0 ) 
+ 4 Re(m~F,Q^- A^_.Q^~,i??*L^)M5SFF5(g., Qt, g, Q, <P~) 
-4mQRe(A^_,Q~,(L?*L^ y,Q<^- + i??*/?^ >^Q,-^-.))M55^^5(g., Qt, g, Q, 0") 
+ AmqRe{m^\*^_~~,\^_Q~MlLT)Y-pssssiq^qs^quQuA~) 

- 4mgRe(m~A^_^-,~~*i?f L?)Y;p555(g, g„ g^, 0") 
+ 2 Re(A^_^_,~~.)Yir555(g, g„ g„ 0") 

+ 4mJmQRe(m~ygQ^-FQg-<^-.i?^Lf )Vc;p^5(g„ g, g, Q, 0") 

- 2 ml{\YQ,^-,Ll\'' + |F^^^_i?^|2)V<.^^^(g„ g, g, Q, 0') 

- 4mgmQRe(FQg<^-*FgQ<^-)V5j,;pP5(gs, g, g, Q, 0") 

+ 4 m,(|y,Q^- 12 + IFq^-^^, nRe(i?^Lf mi)V5^;p^5(g„ g, g, Q, 0") 
+ 4mQRe(m~FQg<^-*FgQ<^-i?f L^)Vcj^j,;P5(g„ ?' ^' '^^ 0") 

- 2 {\Yq,^-*RI\^ + |F,Q^-L^nv5FFF5(gs, g, g, Q, 0")) (13) 



-Re(Ag^g* xx*YFsss(g, g^, g., ^)) 
n?!? = ^-^?fe,. + ^9fc + ^s^iQ„ (14) 

n~|~o = gl{ARe{RlLlY^~,j{Y^,^.~jMsFFSF{qs.q,q,qut^) 

- 4 m,Re(m~y,~o5-.y^o,~ L^L? + mrY~,~{Y^^,.~Rim)MsFFSF{qs. g, g, g~*, Xl) 
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- 2 m~oReim~Y^~oj,Y^~oj,L'iLt + miF^o - l-o,-- i^^i^DM^^^^^lg,, q, q, qt, Xn] 
+ 4 mlRe{Yg~Og*Y^og~Lm)MsppsFiQs, q, q, qt, Xn) 

+ 4m~omgRe(F~o- r~o -Lf i^f + Yg~oj,Yg~Og*LlRl)Mgppsp{qs, q, q, qt, x°) 
-2m~o^mlRe{m^{Y~o^g~Y~o^g~LfLr + Yg~oj,Yg~oj,RiRt)) 

^SFFSFi^s,q,q,qt,Xn) 

-2Re(|y^~o~.|^ + \y^^gqf)yFSSFF{q,qs,qs,q,xl) 

+ 4 m,Re(mi {Y~o,^Y~o,l^RtLr + nk?l^.x05?^ri?f ))Vp,^^^(g, g., gt, g, X° ) 

- 4 m~omqRe(Fg~o~jF~o - )V^s5p^(g, g„ g„ q, Xn) 

+ 4 m~omjRe(F~o - Fg~o~* {m-LlRf + m^LTRl))VpssFF{(l^ Qs, qt, q, Xn) 

- 2 (l^gx°??^^l^ + \^X°nqqt^'^s\'^)^SFFFs{qs, q, q, Xn, Qt) 

+ 4mg(|i;~o~, |2 + |r~o- nRe(m~L^i?f )V5p;p^5(g„ q, q, x^, qt) 
+ 4 rn~oRe{m~Y^~o^~,Y~o^-~LlRl*)Vsppps{qs, g, g, x° , qt) 

- 2 m5(|i;~o5*i?^P + \Y^igg,Ll\'^)ysFFFs(^s, q, q, Xn, gt) 

- 4 m~o mgRe(yg~o~*r~055-J V5^;p^s(g„ g, g, x° , g*) 

+ 4 m~om2Re(miF ~o~.y~o„-~ /2fLf )V5PP^s(g,, g, g, x°, g*; 



(15) 



n~|~_ = ^2(4Re(i?QL?y~_Q.F^_^~)M5FFSF(g.,g,g,g*,Xc") 

^ ^ 9 qxcQl XcQqt '^ t g QXc qt Xc qQs ^ ^' 

^SFFSFiQs,(l,Q,Qt,X7) 

-2m~-Re{m~Y~-^,Ya~-*~*LfLl + mXY*_,~ Y^-^~RfRl) 

Xc ^ 9 qXcQt QXc qt s t g ^^*qQ^ Xc Qqt ^ ^' 

^SFFSFiQs,(l,Q,Qt,X7) 

+ 4 m,mQRe{YQ~-,~.Y^^,_~L',Rf)MsppsFiQs, Q, Q, Qt, Z) 

+ 4m~-mQRe(y~-._^ Y~-a~LfRT + Y ~-^,Ya~-*~,L'iR^) 
Xc ^ ^ Xc qQs Xc Qqt s t q^^ q* qxc qt ^ ^ ■' 

^SFFSFiQs,Q,Q,Qt,Xc) 

-2m~-m„mQRe(mi(F~-*_^ Y--^- 1^*1^ + Y ~-p^,Ya~-*~,RfRl)) 
Xc y V ^ g^ Xc qQs XcQqt « >■ qxcQl QXc qt ^ '- ' ' 

^SFFSFiQs,(l,Q,Qt,Xc) 

~ 2 i\YQ^-*q*f + \YxcQqs\'^)'^PSSFFiq, qs, qs, Q, X7) 

+ ^m,Reim^iY^.^~Y~-Q~RtLr + l^QjrS^Qxr^^^ri^D) 

^FSSFFiQ^9s,qt,Q,X7) 

- 4 m~-mQRe{YQ~-,~.Y~-Q~jVpssFFiQ^ Qs, qs, Q, Z) 
+ 4 m~-m,mQRe{Y~-^~YQ~-,~, {m-LlRT + m^Lf R^) 

'^FSSFFil^Qs,qt,Q,Xc) 
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2 (|F~-5.^'I' + \Y~,,,^m')YsFFFsisis. q, g, Xc", Qt 



- 4 m~-m,Re(r ~-Q.r~-._Q jV5^;p^5(g„ g, g, Xc", Qt) 

+ 4 m~-m2Re(miF ~-^.r~_,_Q^i?^Lr ) V5:p^5(g„ g, g, Xe", Qt)' 



(16) 



A. 2 Neutralinos 

Similarly we find for the neutralinos 

^^'""^^^ " (167r2)2 (%°l?° + %°l?- + ^x°l?+ %«lx« + %?lx-) (1'^) 

where 

%°I0° = -^«^c(4 {-y'^'^m~omlYgg^oRe{X^Og~,Y~OggYg~o~.^)Mg^ 

+ 4 m~oFgg-0oRe(A^o5-^~*F~(,g~Fg~o~.)M55j.i.5(g„ g^, g, g, 0°) 

^ssFFsi^s,quq,qA^) 

+ 4 '^g"^x«^^(^0o?.g;;^WtgS%95t^9?°?l )'^FSSSsiQ^ ^s, qt, qu, (fP) 

+ 2 Re((r^5e°S^5*x% + ^x?^?*^^^-?.) V?=5s V?.c)VF5S5s(g, gs, gt, g«, 0°) 



- 2 mgm~oRe(A^o^o,5-,5*i;^o~.y~o- )Y;p<.g5(g, g„ g^, 0°) 

- ^^{\%0'q,qi (^9X?5j^gx°5r ^ ^x^gg. ^x??? J ) ^^^^^ (^' ^^' ^*' '^°) 
+ 4 (-)''^°'m3m~oi;|<^oRe(Fg~o~.r~o- )V5p^5(g„ g, g, g, 0°) 

+ 2mlY;g^o{\Yg^g,f + |r~o,~ n(V^^^5 + 2(-)l*"IV5^^<,)(g„ g, g, g, 0°) 

+ 4mgm~oFg5<^oRe(yg~o~*F~o- )(2V5^;pps + {-)^^"^^ sfffs){^s, g, g, g, 0°) 

+ 2>;|^c(|i;~o5,p + |F~o,~nV5i.^P5(gs,g,g,g,0°)) (18 

nx°l0- = -^sV(4m,mQm~oRe(FQ5^-,A*_~ ~ rg~o~.FjoQQJM55;p^^ Qi, g, Q, 0") 
+ 4m~oRe(y,Q^-A^_,Q^~.F^OQQ/;x05-.)M55FF5(g., Qt, g, Q, 0~) 
+ 4mQRe(A^_.Q^~,(F~o- F^oqq/q5</'- + ^Qx?Q?^«^<?Q<^-)) 

M55FF5fe,Qt,g,Q,0~) 
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+ 4 m5mQm~oRe(FgQ<^-FQg-^-*F~o- Fg~o~*) V5pppg(g„ g, g, Q, 0") 

+ 4 mgm~o (|FgQ<^- ^ + |>q#-- nRe(V;jo~jF~o- )V5^;p^5(g„ g, g, Q, 0") 
+ 4 mQm~oRe(yQg-^-*ygQ0-yg~o~j>~o- )V5^^;p5(g,, g, g, Q, 0") 

+ 2(|FQ,-^-.y~o-p + |F,Q^-F,~o~.nv5i.i.i.sfe,g,g,Q,0")) (19) 

+Re((y^~o~,y ~o~* + F~o-l~o- )A~~.^^.Yps55(g, g„ g„ X)) 
n~0|~ = n,7r~o|~_^ + 7r~o|~ + 7r~o|Q^ (20) 



+ 4 '^'Re(F^^o~jFJo,-5;F;~o~jr~o- )M55^5^(g,, g, g, g^, x°) 

Msff5f(^., g, ?, ^t, Xn) 

+ 2"^x«"^x?"^9^^((%g-?=^xSlg-5t^gx?5*^gx?5T + ^qxlqt\xlqt^x°qqs'^x'lmt)) 

^SFFsA^s, q, q, qt, xl) 
+ 2 Re((F~o - %g-5 + '^qx°qtKx&t'>^'^qx°q*s^qx°qt + ^x°g-5;^5?°9-5;)) 

VFS5FF(g,gs,gi,g,Xn) 

^FssFFiq^qs,qt,q,x'i) 

^FssFFiq^qs,qt,q,Xn) 
+ 4'^x«^H°Re(^x«5?.^«si?i(^;x?5i^&5; + >;x??r^x«5?J) 

^FssFFiq^qs,qt,q,x'i) 
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+ 4m~om~oRe(y^~o~.y~o -F;~o~jr~o- g, g, Xn, gt) 

+ 4m~omg(|F^~o~,|2 + l>'x«??J')^e(^<?x°5?^x«5?JV5FFFsfe' ^' ?' X°, ?*) 
+ 4m~om~om^Re(Fg~o~*F~Og-~F~o- r^^^^^^ g, g, x°, g*; 
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Xc Xi ^ gxc QJ vxc qt QXiQl 'i^iit Xc qQs XcQqt XiQQs Xiqqt' 

^jFFSpiQs^q^Q^qt^Xc) 

+ 4 m,mQRe(FQ~-*~.K!.._~ F;jo~.r~OQQ JM5^5^(4, g, Q, g„ x,") 

+ 4:m~-mQRe{Y~-,.^Y~-Q~YQ~0Q,Y^^g~ + %-Q.^Qxc-*5?^g*x«g?^x°OQs) 

+ 2 "^x.-"^x«"^'^'^QRe(^xr.-Q7x.-Q5.^Qx?Qjnx?5-r + ^.Xc-Qj%?r?^^x?w7x?55-J 

+ 2 Re((y ~o~ji;~o~. + ^^0-y^o- )(r~-Q~ y~-Q~ + >"Q^-*g*>"Qxc~*9*)) 

VFS5FF(g, gs,gi,Q,Xn) 

+ 4 mom~oRe(y~-^~ Ir-n^r ^c^0;7;r Kc^o;r* + ^r^~-*~.^r)^-*>K^o~.y~o^~ ) 
^ Xi ^ Xc Q9s Xc W<jt Xi 99b 9Xi qt Qxc qt ^Xc ^t qXi qs Xi qqt > 

^FssFFiq^qs,qt,Q,Xn) 
+ 4m~-mgRe(yQ~-*~.r~-^~ (y;o-~r~o,~ + y;^~.y^~o~.)) 

^FssFFiq^qs^qt,Q,Xn) 

+ 4m~-m~om,mQRe(F~,Q~FQ~-.~.(F;~o~.F^o,~ + Y^~,Y^^J) 

^FSSFF(jl^^s,qt,Q,Xn) 

+ 2 (l^9Xc-Q.*^9x^9j I' + l^xc-*5Q/x?.-9J')V5FF^5(gs, g, g, Xc", QO 

+ 4mgm~o(|F~-Q,|' + \Y~^,.Qf)Re{Y^~o~.Y~og~jVsFFFsiQs,q,q,Xc,Qt) 

+ 4^Xc-^x«Re(^«c-Q:^Xc-*5Q/«^$5?JV5i^^F5(^- g, g, x,", 4) 

+ 2^^(inx-Qr^»l' + I^Xc-*.-Q/;x«?/)V5WF5fe,9,?,x;,4) 

+ 4m~-mg(|y^~o~.|2 + \Y~og~f)ReiY^~-Q.Y~-..QjVsFFFsiQs,q,q,x:.Qt) 

+ 4^Xc-^x?^'^e(^9Xc-Q.*^Xc-*5Q/x»5?.^9x??|)V5FFF5(^- ?, g, Xc", Qt 
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(22) 

A. 3 Charginos 

Finally, for the charginos we get 

^^~^ = p^(%ri?° + %n?-+%i5-+%ixo + %rix-) (23) 

with 

+ 2 m~-F,,-^oRe(A^OQ^Q*l~-.5Q/,~-QjM55FFs(Q., Qt, q, q, 0°) 

+ 2m„y„„-^oRe(A,o^ ^JY ~^p,,Y*~~, + (-)\'''"\Y*^,_~ y~-._^ )) 
g gg<t> \ 4>^^QsQ*t^ qXiQt gXi Qt ^ ' Xi gQs Xi gQt" 

^ssFFsiQs^Qt,q,qA'^) 
+ R-((nx-Q.*nkgr + %*.-Q.^r*.-Q>^««^5.^;o.5..5,)Vi^— (^' ^^' ^5. Q., 0°) 

- m,m~-Re(A^o^o.Q^Q.>;5^-Q,>~-*,-Q,)Y:ps5s(g, Q„ Q*, 0°) 

- ^R-(\o,o.5.5.(i;~,5,i;~,5. + ri.^_^7~p^-5j)Y,,,,(g, Q., g. /) 

+ 2 (-)l<^V^m~-F^|^oRe(y ~_^.y~_,_QjV^PP^5(g., g, q, q, 0°) 

+ "^^nio(|F~-Q/ + |y~p,Qj^)(V^^^^ + 2(-)l^"IV^^^^)(a, g, ^^ 

+ 2m,m~-Y^%oReiY^~^^,Y~^,.^J{2VsFFFS + i-)^'^'^'^SFFFs)iQs, ?, ?, 9, /) 
+ n#o(|y~-5/ + |%._^j2)V5i.^p5(4,g,g,g,0°)) (24) 



+ Re(A, ,^ ~.A* ~ ~ (y~-;^.F*~ ~. +y~*_, ~ r~-*_;^ ))yFssssiq,Qs,Qt,qt,(f)~) 

- 2 m,m~-Re(A^_^_.Q^Q.F~-Q,F~-._^^)Y^555(g, g„ g*, 0") 

- R-( V.-*Q.Q, (nk^r^x-Ql + %*,-Q.^--,5;)Y— (^' ^^' ^5. 0-) 
+ 2 m^mQm~-Re(FgQ<^-FQ5^-*r~-,_^^r ~-Q jV5;p^g(g,, g, g, g, 0") 

+ 2 mgmQRe(FQg-^-.FgQ^- )(|y ~_^, P + \y~-*qQ^)'^ sfTps^Q^^ 9' ^' '^^ 0~) 

+ 2 m^m-- (|i;q<^- I^ + \YQq^~, ^)Kft{j^~-^j~-..^)N sp-PFs^Qs. g, g, g, 0") 

+ 2 mQm~-Re(yQg-<^-,ygQ^-F*~_ ~,F^_,_~J V5^^;p5(g„ g, g, g, 0") 

+ (|iQ,-^-*y~-,_Qj2 + |i;o<^-i;~-Q,nv5i.i.i.5(4,g,g,g,0-)) (25) 
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'^x-\x = dnc{2mQm~-Re{Y~-Q~YQ~-,~.^\~~,xx*)^FSssiQ^9s,qt,X) 

+ 2mQm~-Re(rg~0q*F~0QQ,^Qxr*9j-^g*xrQ: ^ ■^QJSlQ?-^x09?.%Q9s%*g-Qt) 

M5FF5F(^.,9,Q,Qt,X°) 
^SFFSFiOs^ Q^ <5, Qt, X°) 

Xn "^ ^ XnQQt XnQQs QXi Qs Xi QQt QXnQt IXnQs qx^ Ql Xi QQs-' 

^SFFSFiQs,q,Q,Qt,xl) 

+ m~om~-m,mQRe((F~o-F~OQQ7Q^-*5-.>;~~Q, + Y^~o^~.Yq~o^.Y~-q~Y~...^J) 

^SFFSpiQs, q,Q,Qt, xi) 
+ Re{(Y~o^-~Y~o^g~ + yqxOj,Y*~oj^)(Y*~-.~YQ~^,~, + Yx; Qq^^^- Qq^)) 

yFssFFiq,qs,qt,q,Xn) 

+ 2mQm~-Re{{Y^o^.~Y~o^g~Y~-Q~YQ~^,~. + Y*~,^~,Yg~o^~.YQ~-,~.Y~-Q~)) 

^TssFFiQ^qs^qt^q^Xn) 

+ 2m~om,Re(yg~Og.F~0 5g,(yQ~-*~.FQ~-.~* + %Q9,>~-q5)) 

^ Fss-FF^Q^^s^quq^xl) 
+ 2m~om,mQm~~Re(y~o,-5;y,jo^.(%P~jl^^^ 

^Fss-FFiQ^qs,qt,q,Xn) 

+ 2m,m~.(|F^~o5*P + \Y~o^^~f)ReiY^~^Q^Y~^,.^jVsFFFsiQs.q.q.xlqt) 

+ 2 '^x»^xr^^(^</x^5-r%5?*^;xrQj^^-*5Q J^^^^^^^^- ^' ^' ^° ' ^*) 

+ ^'(l^5X$i5?%-5Qj' + \'^X?.^qXqxTQ'f^^SFFFsiQs^ 5' ?' X" , ^t) 

+ 2m~om,(|r~_^.|2 + |l~_._^j2)Re(F^~oj.F~o-)V5^^5(Q,,g,g,X°,gt) 

+ 2 "^x»^x-^'^e(^5x°?^x°^?t^x-*g-Q.%-Qp^SFFFs(<3., g, g, X° , qt)) 

(27) 

x^\x ^s cy \\ g^^Q, q^^Q* X, qQt X, qQs'^ Xc gQs Xc qQt qxcQt iXaQl" 
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VFSSFFiq,Qs,Qt,q,Xn) 

^ Xi ^ Xc gQs Xc gQt Xi qQs QXi Qt iXcQt 'iXcQt QXiQt Xi iQt' 

^FSSFFiQ^Qs,Qt,q,Xn) 

Xc 1 ^ qXcQt Xc gQs^ Xi gQs Xi gQt gXiQt iXiQt" 

^ FSSFFy^i QsiQti Q, Xn ) 

+ 2m~-m~-m„m„Re(Y~-,.p,Y~-p^jY*~^~Y*^,_~ +Y ~-x,Y~-,.p, )) 
Xc Xi 1 1 ^ Xc gQt gxcQi^ gxiQt Xi gQt gXi Qt Xi gQs" 

^TSSFF\^i QsiQti Q, Xn ) 

+ ^KxcQ*/gx-Q*s\' + l^xr.-o/xr*.-QJ')^^^^^^^^- ^' ^' ^^"' ^'] 

+ 2 m,m~- (IF ~_^. p + |F~_..^J^)Re(r ~-^.r~-._^jV,^^^,(Q., g, g, x.", Qt 

+ 2 -,~,-^-Re(y ~,5.y~^_.^_5 Fl_ ~.Kr..^.~ jV,,,^,(Q., g, g, x^, QO 

+ ^S(lnx.-Q:%*.-Qj' + \'^x-gQt'^gxTQt\'^'^SFFFsiQs,<l,Q,X::.Qt) 

+ 2 "^Xe-^.d^x-Q^ I' + \^xr^Qs\")^<^gx,Q*/x-gQt)'^SFFFsiQs, <1, 1, Xj, Qt 

+ 2"^x.-^x-^?R^(^.xc-Qr^xr.-Q/x-*.-Q7.x-gP^^^^^(^- ^' ^' ^^"' ^*)- 



A. 4 Yukawa couplings 

The relevant couplings are 



YuuhO = 


-1=008 a Yu 

v2 


YnuAO = 


COS pYu 


V - — 


: ^sm(3Yd 


^ddhO = 


- V2''''''^' 


^ddAO = 


^sm/3r. 


YduH~* = 


cospYu 


yuuHO = 


i=smaYu 

V2 


yuuGO = 


sm(3Yu 


^udG" = 


: ^cosPYd 


YddHO = 


-^ cos aYd 

v2 


YddGO = 


1= cos p Yd 

v2 


YduG— = 


sm f3Yu 



(29) 



and 
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A 



h°h°usu* 



A 



hOhOdsd* 
HOROZsUt 



HOHOdsd* 



A 



H~ H~*Usut 



-6s,t sin^ a Y^ 
-6s,t sin^ a Y^ 
-6s,t cos^ a Y^ 



AOAOusU* 



AOA^dsd* 



'GOGoZd* 



G-G—dsd* 



sin^ /? L^LTY] - cos^ /3 R^RTY^ 
sm^ P RiRtY! - cos-" (3 LiLfY^ 
cos' (3 L^LrY,'- sin' (3 R:RrY^ 
cos^ /3 RiRf*Y,' - sin^ /? Lf Lf i;^ 



-5s,tCos'pY^ 
-5s,tsm''pYj 
-Ss,tsm'(3Y;^ 
-6s,t oos' (3 Yl 



(30) 



The four squark couplings are 



UsU-UuUZ 



UsU^dudl 
dsd^UuUl 
dsd'dud* 

dsd*du,gd*g 



-{R^RrLm* + l:ltrirt)y^ 

T d T d* -QU ■nu*\/-2 zjd rid* tut u*\^2 

— 1j„1j+ -ft,, -ft,, -f,, — rL„rLf 1j,,1j,, I j 



s t -^^u ^v u 



s H u V 



A- 



-{RiRf^LtLt + LiLt'-RtRtnY,' 

-{R:LrLi^Rr + L:RrRi^Lr)YuYu, 

A~~~ ~. =0 
-{RiLfL'tRt^' 



LtRTRi'LtnYdY. 



and finally the Higgs couple to squarks according to 



A 

A 
A 

A. 



H~ Usdt 



m*Limu cos (3 + RrHAu cos (3 + jj,* sin (3)1^ + rudRi cos /3))K 
-{R^Rimu sin p + L^ (i?f (// cos p + A*^ sin /3) + m^Lf sin /3))rd 



A 



// utd* 



G~ Usd* 



(31) 



m*Lim^ sin /3 + i?r((^« sin (3 - jj.* cos /5)Lf + m^i^f sin (3))Y^ 
-{RTRimuCos^ + Lr(i?f(-/isin/3 + A^cos/3) + m^Lf cos /5))lrf 



A* ~ ~ 

G^utd* 



h'^UsU* 



2mu6s,t cos a + L^R^* {A^ cos a + //* sin a) 

+L'l*R'^{fi sin a + a; cos a)) r„ 



A 



h°dsd* 



V2 



—j=\2nid5s,t sin a + L'^R'l:*{Ad sin a + //* cos a) 

V 2 ^ 



A 



HOusU* 



+L'l*R'l{ii COS a + A^ sin a)jYd 
-j= \2mu5s,t sin a + L^Rt* {Au sin a — n* cos a) 



v^ 
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+LJ'*i?^(A* sin a - /i cos a) j F^ 
A^ojj* = -^ ( 27nd6 s, t cos a + L'^Rf*{ Ad cos a — fi* sin a) 



+Lf i?f (A^ cos a - /i sin a) jF^ 
A^o^^~, = -^{l: Rr i A^ cos (3 + fi* sin (3) -LrR:iA: cos P + fi sin P))y^ 

X^odJt = -^{LiRf*iAdsinP + fi* cosf3) - Lf*Ri{A:sinP + ficosp)Yd 

Ago^^~. = -^{L^RriAuSinP - i^* cosP) + LrR:{fxcosP - A:sinP))Y^ 

\goIJ* = ^ {L^Rf* (/i* sin /5 - A, cos /3) + Lf /2f ( A* cos /5 - /i sin /?) F, 

(32) 

In Eq. ( I29]|32l) u and rf denote the up-type resp. down-type quark, the generation 
index is suppressed. 
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Figure 6: Left: One-loop (red) and two-loop (solid) neutralino pole masses as a function 
0/ tan/3. Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQCD 
corrections only, the solid lines also the Yukawa corrections. 
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Figure 7: Left: One-loop (red) and two-loop (solid) chargino pole masses as a function 
o/tan/?. Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQGD 
corrections only, the solid lines also the Yukawa corrections. 
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Figure 8: Lejt: One-loop (red) and two-loop (solid) neutralino pole masses as a function 
of At. Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQCD 
corrections only, the solid lines also the Yukawa corrections. 
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Figure 9: Left: One-loop (red) and two-loop (solid) chargino pole masses as a function 
of At. Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQGD 
corrections only, the solid lines also the Yukawa corrections. 
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Figure 10: Left: One-loop (red) and two-loop (solid) neutralino pole masses as a function 
of (pAf Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQCD 
corrections only, the solid lines also the Yukawa corrections. 
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Figure 11: Left: One-loop (red) and two-loop (solid) gluino pole mass as a function o/0m3- 
Right: Two-loop mass-shift in GeV. The dotted lines contain the SUSYQCD corrections 
only, the solid lines also the Yukawa corrections. 
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Figure 12: Left: two-loop pole mass of Xi os a function of (pM^- Right: Xi ■ The dotted 
lines contain the SUSYQCD corrections only, the solid lines also the Yukawa corrections. 
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